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Abstract 

A two-sex Basic Reproduction Number (BRN) is used to investigate the condi- 
tions under which the Human Immunodeficiency Virus (HIV) may spread through 
heterosexual contacts in Sub-Saharan Africa. (The BRN is the expected number 
of new infections generated by one infected individual; the disease spreads if the 
BRN is larger than 1). A simple analytical expression for the BRN is derived on 
the basis of recent data on survival rates, transmission probabilities, and levels of 
sexual activity. Baseline results show that in the population at large (characterized 
by equal numbers of men and women) the BRN is larger than 1 if every year each 
person has 82 sexual contacts with different partners, the BRN is also larger than 
1 for commercial sex workers (CSWs) and their clients (two populations of different 
sizes) if each CSW has about 256 clients per year and each client visits one CSW 
every two weeks. A sensitivity analysis explores the effect on the BRN of a dou- 
bling (or a halving) of the transmission probabilities. Implications and extensions 
are discussed. 

KEYWORDS: Basic reproduction number, transmission probability, log-log com- 
plementary model, Weibull distribtution. 
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1 Introduction 



There is a growing debate as to whether hete rosexual contacts are the main mode of 
transmission of HIV-l in Sub-Saharan Africa ( 



2003 



Brewer et al. 



2003 



Gisse. 



French. Riley and Garnett 



auist. 



2006 



'otterat. Brodv and Vachon 



Deuchert and Brody 



2007h . 



Some question the conventional wisdom of a heterosexual epidemic on the basis tha t 



20071). 



assumed transmission probabilities per coital act are inflated (iDeuchert and Brody 
Others, on the contrary, have developed complex mathematical models to show that 
alternatives such as the use of unsafe medical injections were unlikely to be the main route 
of transmission bec ause they rely on unfeasibly high iatrogenic transmission probabilities 



( iFrench et al 



20061) 



In this paper we will show that a complex mathematical model is not necessary in order 
to assess the feasibility of a heterosexual epidemic. Indeed, the question can be studied 
by focusing on the basic reproduction number (BRN), which is the expected number 
of secondary infections generated by one infected individual in a completely susceptible 
population (i.e. at the beginning of an epidemic). The disease will spread if and only if 
the BRN is > 1, i.e. each infected individual infects more than one other person. 

The calculation of the BRN hinges crucially on the evolution over the course of the 
infection of the transmission probability per coital act. A careful study of a population- 
based cohort of discor dant couples (one person infected) i n Rakai, Uganda, has shed 



light on this question (Gray et al. 



2001 



Wawer et al. 



20051 ). The first high- infect ivity 



stage of the infection is characterized by an early peak in the viral load. This pattern is 
paralleled by a rise in the transmission probability per coital act that reaches a peak of 
about 0.008 before declining sharply one year into the infection. During the long second 
(asymptomatic) stage, the viral load is very low. The probability of transmission remains 
also very low, at around 0.001 per coital act. The third and last stage of the infection is 
characterized by a late peak in the viral load (and in the probability of transmission). 



The likelihood of a heterosexual epidemic depends on the transmission probabilities 
but also on the number of partners. For example, with such relatively low probabilities, 
the disease may not take hold in a serially monogamous population, but could spread 
within high activity groups characterized by rapid changes in partnerships. 

For this reason it is sufficient to focus on the possible heterosexual spread between 
high-activity groups. Indeed, if the disease can spread between such groups, it will spill 
over to others even in the transmission is inefficient from high to low activity sexual 
partners. The example that comes to mind is that of an epidemic that may spread 
efficiently between commercial sex workers (CSWs) and their clients. The latter, in turn, 
may infect, however inefficiently, their long-term female partners and thus spread the 
virus significantly among low activity women. 

We will first derive an analytical expression for the viral load on the basis of re 



cent 



y available information obtained from the Rakai study (iGray et al. 



2001 



Wawer et al. 



20051 ). We will then use the log-log complementary model to obtain an expression for the 
probability of transmission per coital act. This probability will be a function (via the 
viral load) of the infective age ia (time since infection) and of the infective age at death 
iad (time from infection to death). We will call ptr{ia,iad) this transmission probability. 

Highly active groups (such as CSWs and their clients) will be characterized by annu- 
alized numbers of coital acts NCA{ia,iad) (assumed to take place always with new part- 
ners). This number depends on the infective age and the infective age at death because 



the number of coita 



death fi Wawer et al. 



acts decreases as a person advances in the disease and approaches 



20051). 

The transmission rate for an individual who has been infected ia years and will die at 
iad years is now NCA{ia,iad)ptr{ia,iad). If s{x) is the density function of the infective 
age at death, the basic reproduction number Rq is the expected value of the number of 
secondary infections generated by one individual during his/her infective life course (with 



3 



maximum duration u): 



def. 



s{y) / NCA{x,y)ptr{x,y)dxdy. 



This is the BRN of a single sex model in which transmission is the same between 
all individuals. In particular the transmission rate and density of survival time are as- 
sumed to be the same for b oth sexes. However, survival differs slightly between men 



and women (UN AIDS 



2002). There can also be large differences in the number of 
coital acts when dealing with groups of different sizes such as CSWs and male clients. 
Transmission rates may be different for the two sexes. For example there is a grow- 
i ng consensus that male circumcision red uces significantly female-to-male transmission 



( iNagelkerke. Moses, de Vlas. Bailey 



20071 ). 



For these reasons we must consider two sex-specific basic reproduction numbers: the 
number Rf^ of secondary males infected by one infected woman ("the female to male 
BRN"), as well as Rmf, the male to female BRN. The product RfmRmf is therefore the 
number of same-sex tertiary infections generated by one infected person. The expressions 
for Rmf and Rfm will be those of Eq. ([T]) with the three functions s{y), NCA{x,y) and 
ptr{x,y) indexed by m for Rmf and by / for Rfm- 

Although one could have defined the two-sex BRN as the product RfmRmf, it is 
generally defined as the harmonic mean 



def. 



\/ RfniR 



mf 



(2) 



of the two sex-specific BRNs. This definition re: 



over two generations (IHeesterbeek and Roberts 



ects the fact that transmission takes place 



20071 ). It is also consistent with the defini 

-Km/ 



tion of the BRN as the dominant eigenvalue of the next-generation matrix 





Rfm 
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The threshold condition for an epidemic flare-up is now Rq > 1, i.e. RfmRmf > 1: 
the number of same-sex tertiary infections generated by one infected individual must be 
larger than 1 for the epidemic to take hold. 

We will show that with a set of realistic baseline parameter values, then in the popula- 
tion at large (characterized by equal numbers of men and women) , the basic reproduction 
number Rq is larger than 1 if each year every person has 82 sexual contacts with different 
partners. Within the CSW-client populations (where men outnumber women) the infec- 
tion can spread if each CSW has about 256 clients per year and each client visits one 
CSW every two weeks. 

The paper is organized as follows. In Section 2 below we derive expressions for the 
viral load, the transmission probability per coital act, the annualized number of coital 
acts, and the density function of survival times. In Section 3 we give an expression for 
the two-sex basic reproduction number and formulate the threshold condition in terms of 
the Index of Sexual Activity. Results are then illustrated with realistic parameter values 
pertaining to Sub-Saharan Africa. The sensitivity of the results are discussed for different 
values of the probability of transmission function. In Section 4 we discuss our findings, 
their implications and possible extensions. 

2 The four components of the basic reproduction 
number 

The construction of the basic reproduction number is the same for both sexes. For ease of 
exposition we will therefore drop the indexes / and m from the functions (and parameters) 
used to define Rfm and Rmf- 
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2.1 Viral load 



In the absence of treatment, the logarithm base 10 (Ioqiq) of the viral load (measured in 
copies/mL) follows a well-establi shed pattern as a function of the infective age ia and of 



the infective age at death (iad) (IRapatski. Suppe. Yorke 



20051 ). During the first year of 



infection the logarithm increases to approximately 5 and decreases rapidly thereafter (first 
stage). It then remains around 3 during the long asymptomatic second stage. About a 
year before death there is a second peak in the viral load. 

We now describe a function noted LVl{ia, iad), that captures this "twin peaks" pattern 
in the Ioqiq of the viral load during the course of the infection. The parameters that define 



the function are given in Table 1, together with b aseline nunierica. 



empirical results obtained from the Rakai study (IGray et al. 
We take the same parameter values for both sexes. 



2001 



values which reflect 



Wawer et al. 



20051). 



Table 1: Parameters of LVl{ia,iad) 


Parameter 


Baseline value 


iai'. Infect, age at first peak 


0.4 year 


Mi: Value of LVl{iai,iad) at 1st peak 


5 


m: Low value of LVl{ia,iad) during second stage 


3 


Ti'. Time preceding death at 2nd peak 


1 year 


M2: Value of LVl{iai,iad) at 2nd peak 


4.8 


ai,a2'. Parameters that determine the variance in the 1st peak 


1.3; 0.2 


03: Parameter that determines the variance in the 2nd peak 


0.7 



We first define the function 



hi{x, ai, Ml, iai 



def. 



Mix"^ exp 



x(l — ai) 



lai 



—1 



X exp{l — ai] 



(3) 



which reaches a maximum of Mi for x = iai. We then need the largest root x* of the 
equation hi{x, ai, Mi, iai) = rn in the unknown x; x* = 1.647 and will be used to obtain 
the low value m during the long asymptomatic second stage. 
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We also need the function 



def 

h2{ia,a2,x*) =' + exp(— 0:2)] x 



1 + exp a2 



ia{l + exp{a2)) 



'1 + exp{a2)) ^ 



(4) 



which will be used as the argument x in the function hi. This will produce the first peak 
followed by the low value during the asymptomatic stage. 
The function 

def. 



/i3(x, y, ^3, Ti) =■ exp [-asix - y + Tif] 



(5) 



will be used to obtain the late-stage peak. We combine these elements to finally define 



def 

LVl{ia, iad) = hi{h2{ia, 02, a;*), ai, Mi, iai) + 

[M2 — hi{h2{ia, a2, x*), ai, Mi,iai)]h3{ia, iad, as, ri). (6) 



The loQiQ of the viral load function LVl{ia, iad) corresponding to the parameter values 
in Table 1 is plotted in Figure 1 for an infective age at death (iad) of 7 years (together 
with the transmission probability function derived below). For a later infective age at 
death the function is similar with just a longer asymptomatic stage. 



2.2 Transmission probability per coital act 



Gray et al. 



(120011 ) propose the log-log complementary model for the transmission prob- 
ability per coital act with (chronological) age and viral load as covariates (variables age 
and vl). Under this model this probability of transmission is of the form 



def 

ptrQ^age, vl) =' 1 — exp[—exp{KQ + kivI + K2age)] 



(7) 
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Figure 1: logio of viral load function LVl{ia, iad) and 1000 times probability of transmis- 
sion per coital act function (1000ptr(ia, mrf)) for a person who dies seven years into the 
infection {iad = 7). 



with parameters K,m{fn = 0, 1, 2). 



The effect of age does not appear to be very strong (iGray et al. 



2001 



Wawer et al. 



20051 ) and would complicate the expression for the basic reproduction number. As a 
simplification we therefore drop age as a covariate. Bearing in mind that iQLVi(ta,tad) 
the viral load, we re-express an average (across ages) probability of transmission per coital 



act as the function 



ptr{ia, iad) 1 - exp[-exp{Ko + k^IO^^^'^'^''^'^^)]. (8) 

We parameterize this function by specifying the values ptrhi and ptrio of ptr{ia,iad) at 
the values Mi and m of LVl{ia,iad) corresponding to the first peak in viral load and to 
the low plateau. For given values of ptr^i and ptrio, the parameters Kq and Ki are then 
the roots of the system 

ptrhi = 1 — exp[—exp{Ko + KiIO^^)], ptrio = 1 — exp[—exp{Ko + kiW^)] (9) 
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from which 



In 



ln{l - ptrio 
ln{l -ptThi 



In 



IQMl-m _ I 



+ In [ln{l - ptrio) ^] , Ki 



ln{l - ptrio) 
ln{l -ptThi) 



IQA/i 



(10) 



With the numerical values ptrio = 0.001 and ptr^i = 0.008, the resulting function 
ptr{ia, iad) (multiplied by 1000 in Figure 1) provides a good stylized ap proximation of 



recent empirical estimates based on the Rakai study (IWawer et al 
same parameter values for both sexes. 



20051 ). We take the 



2.3 Annualized number of coital acts with different partners 

We next construct a functional form for the annualized number of coital acts NCA{ia^ iad). 
This function w ill reflect a decreas ing level of sexual activity as an infected person ap- 



proaches death (IWawer et al. 



20051 ). The parameter A will be the value of NCA at the 
time of infection {NCA{Q^iad) = A), i.e. the annual number in the absence of HIV 
infection. 

The parameter will be the fractional number of coital acts remaining when an 
individual reaches the infective age iad — ti at which the viral load reaches its second 
(pre-death) peak {NCA{iad — Ti,iad) = A0). Finally NCA will be at the time of death 
{NCA{iad,iad) = 0). 

We now define for ia < iad the function G{ia, iad) equal to the fractional number of 
(annualized) coital acts remaining for an individual infected ia years ago and who will die 



iad years into the disease: 



G{ia, iad) 



def. 



1 — ia/iad 



ia{Ti — (p.iad) 
iad.(f){iad — ti) 



if iad > Ti, 



if iad < Ti. 
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Figure 2: Fractional number G{ia,iad) of (annualized) coital acts remaining ia years into 
the infection for an individual who will die iad=l.l, 3 or 5 years into the infection (with 
ri = 1 and = 0.61, the function satisfies G(0.1, 1.1) = G(2, 3) = G(4, 5) = = 0.61). 

This function is equal to 1 for ia = 0, to for ia = iad—Ti and to for ia = iad. A function 
NCA{ia, iad) that has the required properties is obtained by multiplying G{ia, iad) by 
A: 



def. 

NCA{ia,iad) =' A.G{ia,iad). 



(12) 



A baseline value of was taken equal to 0.61 on the basis of a mean reported number 
of coital acts per week of 10.2 at the beginning of the infection and of a mean number 



20051), 10.2/6.2 = 0.61). 



during a 6-15 month period prior to death of 6.2 ( (jWawer et al.l 
We take the same parameter values for both sexes. 

The function G{ia, iad) is plotted in Figure 2 for three different values of iad and with 
Ti = 1. As iad becomes closer to ri, the function G{ia, iad) of ia approaches more and 
more rapidly as ia tends to iad. The fact that G{ia,iad) is zero when the infective age 
at death iad drops below ti means that no sexual activity is assumed for a very short 
infection (e.g. an infection that lasts less than one year when ri = 1). This drop to zero 
in sexual activity may not be entirely realistic, but is of little importance since there are 



extremely few, if any, infected individuals who will survive such a short period. 
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2.4 Density function of infective age at death iad 



Following the Wo rld Health Orga nization we assume a WeibuU distribution for the infec- 



tive age at death flUNAIDS 



2002I ). We parameterize this distribution with its median me 
and shape parameter /3. If we define a me {ln{2))~^^^ the density function s{x) of iad 
is then 



S[X] 



-exp 



x\f^ 
a 



(13) 



We use (3 = 2.5 for the shape parameter for both sexes and a slightly shor ter median 



for women (me/ = 8.6 years) than for men {mem = 9.4 years) (jUNAIDS 



2OO2I). 



3 Results 



3.1 Threshold conditions on the basic reproduction number 

In general all functions (and parameters) are indexed by / and m. The two sex-specific 
basic reproduction numbers are then 



def n 

Rfm ='A/ / Sf{y) / Gf{x,y)ptrf{x,y)dxdy 

Jy=0 Jx=0 



(14) 



Rm.f — A 



y 



(y) / Gmix,y)ptrmix,y)dxdy. 

y=Q Jx=0 



We next define the quantity 



(15) 



Io'^= ( / Sf{y) / Gf{x,y)ptrf{x,y)dxdy X / Sm{y) / Gm{x,y)ptrm{x,y)dxdy 

\Jy=0 Jx=0 Jy=0 Jx=0 

(16) 

This quantity /q reflects at the individual level the combined effects for both sexes of 
variable infectivity, mortality, and sexual activity over the course of the infection. We 
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also define the Index of Sexual Activity (ISA) as the harmonic mean of the contact rates 
Am and Af (i.e. the sex-specific annualized numbers of coital acts at the beginning of 
the infection): 

ISA v/a:::a7. (i?) 

The Index of Sexual Activity measures the level of sexual activity between the two groups. 
With these definitions, the composite basic reproduction number 

Ro — \/ RfmRmf (18) 

will be larger than 1 if and only if 

ISA > lo. (19) 

The (annualized) number of coital acts men have with women must be the same as the 
number of acts women have with men. If Pf and are the sizes of the corresponding 
female and male populations, we must therefore have 

PfxAf = PmxA^. (20) 

The threshold condition (fT9!) can then be paraphrased by saying that when both popula- 
tions have the same size then Jq is the minimum annualized number of coital acts each 
person must have with different partners in order for the disease to take hold (since then 
Af and are equal). 

3.2 Numerical illustration (with baseline parameter values) 

The baseline parameter values and functions are those given in Section 2. They are the 
same for both sexes, except 
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50 100 150 



A 

m 

Figure 3: Phase space of annualized numbers of coital acts by women (A/) and men 
{Am) with locus Am A/ = Jq of values for which Ro =1 (baseline hyperbola obtained with 
parameter values of Section 2). The basic reproduction number Rq is larger than 1 above 
the curve and vice-versa. The black circle on the hyperbola is the fixed point /q = 81.60, 
i.e. the minimum annual number of coital acts for the epidemic to spread when the male 
and female populations are of the same size. The "feasible rectangle" covers a range of 
plausible values of A^ and Aj for commercial sex workers and their clients (see text). 
The basic reproduction number at the four corners show that Rq is larger than one in 
almost the entire feasible rectangle. The two other hyperbolae correspond to a halving 
and to a doubhng of the probabihty of transmission function ptr for both sexes. The 
effect is linear on Jq which is then 40.8 and 163.2. With a halving of the probabilities, Rq 
is larger than one only for high levels of sexual activity (for example Am = 100 prostitute 
visits per year and A/ = 400 customers per year for each commercial sex worker). With 
a doubling of the probabilities, Ro is larger than one in the entire feasible rectangle and 
well below. 
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• for the crucial sex-specific contact rates and Af between the two groups that 
will be used for the sensitivity analysis below. 

• for the slightly different median survival times me/ = 8.6 and mem = 9.4. 

The quantity Iq is independent of and Aj and its baseline value is 81.60. In the 



(Am, A/) phase space the corresponding baseline hyperbola of equation A^Aj 



is 



therefore the locus of values for which _Ro is equal to 1 (Figure 3). 

With male and female populations of the same size (A^ = Aj) the fixed point Iq = 
81.60 (black circle) of the baseline hyperbola tells us that each newly infected person 
needs about 82 coital acts per year with different partners in order for the disease to 
spread. This is a small number compared to the documented 10.2 acts per week reported 
above. The requirement that these acts take place with different partners, on the other 
hand, is in stark contrast w ith surveys that report an average of about one partner per 



year in Sub-Saharan Africa ( jPeuchert and Brody 



20071 ). In short, the disease can spread 



between groups of men and women of equal sizes for a reasonable annual number of coital 
acts, but with the requirement of a very high turnover of partners. 

Commercial sex workers and their clients are groups of different sizes characterized 
by a (Am, Aj) point that lies above the main diagonal (Am=Aj) of the phase-space 
diagram in Figure 3. The number of acts per year varies considerably, however, with 
estimates for men in th e range 0.5-2 prostitute visits per week, i.e. 26 < A^ < 104 



(jNagelkerke et al. 



20071 ). A range for annual numbe rs of clients is based on medians of 



4 and 9 per week in rural and urban areas of Kenya (lElmore-Meegan. Conroy and Agala 



20041 ). These medians translate into the range 208 < Aj < 468 for annualized numbers of 
coital acts performed by each CSW with her clients. The lower bound may reflect " casual" 
practices, while the upp er one is probably conserva tive, with estimates of up to 15 per day 



(Aj = 5475) in Ghana (lAsamoah-Adu et al 



200ll ). The resulting "feasible rectangle" of 



(Am, A/) values is depicted in Figure 3. The values of the basic reproduction number _Ro 
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(Eq. (ITSl) ) at the four comers are given in the rectangle. At the lower left corner each 
infected CSW will infect Rfm = 2.47 clients, who will each in turn infect Rmf = 0.33 



CSWs. The resulting Rq is equal to ^JRfmRmf = 0.90. This shows that despite an 
efficient female-to-male transmission the infection will not spread if clients visit a CSW 
only every other week and CSWs service only 208 clients a year. If (A^, A/) moves up 
the left side of the rectangle then Rmf remains unchanged and the point (A^, A/) enters 
the Rq> I region for A/ = II/ \n = 81.60726 = 256.1 (i.e. when Rfm reaches 3.04). 
For (Am, A/) at the upper right corner, Rfm = 5.55 and Rmf = 1.32: both BRNs are 
larger than 1 for an overall Ro of 2.70. 

3.3 Sensitivity analysis on transmission probabilities 

There is a fair amount of uncertainly concerning the values of the peak and low transmis- 
sion probabilities ptrhi and ptrio whose baseline values were taken as 0.008 and 0.001. For 
example the 0.008 figure was an estimated probability at infe ctive age 5 months, with a 



95% confidence interval of (0.004, 0.0015) fWawer et al.ll2005h . 



In order to assess this sensitivity to the transmission probabilities we plotted in Figure 
3 the two hyperbolae corresponding to a halving (resp. a doubling) for both sexes of both 
parameters ptvhi and ptrio- This means a halving (resp. a doubling) of the ptr function 
and has a linear effect on Jo which is doubled to 163.2 (resp. halved to 40.8). With 
a halving of the ptr function about two thirds of the feasible rectangle is in the Rq < 1 
region. Fairly high contact rates Am and A f are needed in order for the disease to spread. 
With a doubling of ptr the rectangle is entirely in the Rq < 1 region and the spread will 
take place with much smaller values of Am and Af. 
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4 Discussion 



Our goal was to use the basic reproduction number Rq to investigate whether HIV-1 



can spread in Sub-Saharan Africa (or e 



se wh ere) primarily through heterosexual contacts. 



Fraser. Riley. Anderson and FergusonI (120041 ) point out that published estimates of Rq for 



generalized heterosexual HIV epidemics are hard to come by. This is no doubt because 
of complex transmission mechanisms and the heterogeneity of the populations involved. 

In this paper we have introduced a data-driven analytical expression for a two-sex basic 
reproduction number that captures the nuances of variable infectivity and of changing 
levels of sexual activity over the course of the infection. 

The results show that with our baseline parameter values the disease can spread be- 
tween groups of men and women of equal sizes with a reasonable number of coital acts per 
year (82). However these contacts must take place with different partners, which implies 
an unusually high level of sexual activity. The infection will spread between CSWs and 
their clients for most plausible contact rates and Af. 

By focusing on the BRN we were able to investigate conditions for a heterosexually 
driven HIV/AIDS epidemic on the basis of the survival, infectivity, and contact rates 
alone. We did not have to make the many complex assumptions needed for a full-blown 
dynamic model of HIV transmission (e.g. mixing patterns, partnership formation rules, 
etc.). Valuable insights have thus been gained with a minimum number of assumptions 
and parameters. 

Epidemiologists, public health officials and others can now make judgements concern- 
ing the potential spread of HIV-1 by checking the simple condition A^Aj > Jq with 
local values of the contact rates A^ and Af. It is wise however to take into account the 
uncertainty on Iq by considering its plausible lower and higher values of 40.8 and 163.2. 

At least the overall (average) transmission probability obtained from the Rakai study 
is believed to be similar to that reported from "prospective studies of European, north 
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American and Thai heterosexual couples" (IGray et al. 



2001 



p. 1152). The transmission 



function ptr may therefore be applicable to populations outside of Uganda. 

Our results assume no intervention (condoms, circumcision, therapy, etc). However 
the survival, viral load or transmission functions can be changed to assess the effect on 
the basic reproduction number of longer survival and/or antiretroviral therapies. More 
generally, the theoretical results of Section 3.1 can be used for any sexually transmitted 
infection by using the appropriate function needed to calculate Rfm and Rmf- 

Finally we emphasize that our results are based on the assumption of instantaneous 
partner changes which implies a high level of sexual activity. Our results are not appli- 
cable to a population characterized by alternating periods in and out of a partnership 



as m 



Morris and Kretzschmar 



2000l ). Everything else being equal (survival, infectivity, 
etc.) expected numbers of secondary infections will be lower in such partnerships. This is 
because during the whole duration of a (monogamous) partnership an infected individual 
can transmit the virus to only one person. If that happens, there will be no more trans- 
mission until the end of the current partnership, itself often followed by a period without 
sexual activity. Calculating an expected number of secondary (and tertiary) infections in 
this case is a difficult statistical problem currently under investigation. 
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